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A new method to constrain the cosmological equation of state is proposed by using combined samples of gamma- 
ray bursts (GRBs) and supernovae (SNela). The Chevallier-Polarski-Linder parameterization is adopted for the 
equation of state in order to find out a realistic approach to achieve the deceleration/acceleration transition phase 
of dark energy models. As results, we find that GRBs, calibrated by SNela, could be, at least, good distance 
indicators capable of discriminating cosmological models with respect to ACDM at high redshift. 



1. INTRODUCTION 

The acceleration of the universe is one of the 
most important scientific discoveries of the recent 
years, V r ^ led to the hypothesis of the existence 
of a force, called dark energy, as responsible for 
this acceleration. The observational evidence for 
an accelerating universe came from the analysis of 
the luminosity distances of Type la Supernovae, 
obtained by the well-known relationship between 
the luminosity and the light curve shape, [2]. So 
one of the main research lines in observational 
cosmology is to reconstruct the evolution of the 
universe not only from a dynamical point of view, 
but mainly from a thermal one. 

On the other hand the limits imposed by the 
instruments devoted to the detection of Super- 
novae la do not allow us to go beyond a certain 
redshift value, fixed by the redshift of the most 
distant supernova yet seen, z t = 1.7. Then in 
order to measure distances beyond this value of 
redshift, we need indicators observed at higher 
redshift than z t . One of the possible solutions 
to this problem is the use of Gamma Ray Bursts 
(GRB). GRBs are the most powerful explosions 
in the universe, since they originated from the 
formation of a Black Hole. Moreover they are ob- 
served at considerable distances, of the order of 



the redshift z — 8 — 10, so that there are sev- 
eral efforts to frame them into the standard of 
the cosmological distance ladder. 

In the literature there are several detailed mod- 
els which give account for the GRB formation and 
emission mechanisms, but none of them is intrin- 
sically capable of connecting all the GRB observ- 
able quantities, see e.g. [3]. For this reason GRBs 
cannot be considered as standard candles. How- 
ever there exists several observational correlations 
among some photometric and spectroscopic prop- 
erties of GRBs which features them to be used as 
distance indicators. 

In this work we consider two empirical correla- 
tions between GRB's observable quantities in or- 
der to build a GRB Hubble diagram, after a cali- 
bration of these correlations with the Supernovae 
la data. This GRB Hubble diagram represents 
the starting point for our analysis, in which we 
take into account the possibility that GRBs may 
be an extension of the Supernovae la at high red- 
shifts. Later on the analysis of this new sample 
is made in order to test the Dark Energy models 
in the context of the ACDM concordance model. 
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2. THE GRB HUBBLE DIAGRAM 

In recent years, a great number of space mis- 
sions have been devoted to detailed and con- 
tinuous observation of the GRBs. In particu- 
lar, among the many results obtained, there were 
found some relations between spectroscopic and 
photometric quantities of GRBs themselves to the 
point that, using these relations, it is possible to 
trace back a very accurate estimate of the total 
energy emitted by these sources. Nevertheless, up 
to now, there is no theoretical model that fully ex- 
plain these relations so for this reason GRBs can- 
not be considered as standard candles in a proper 
sense. For a detailed review of the observational 
features see @]. 

In this work we are taking into account the ex- 
isting 3-parameter relations. These relations put 
the better constraints on the data giving less scat- 
ter between the theoretical relation and the ex- 
perimental data. The first relation is the Liang- 
Zhang relation, [5], which connects the isotropic 
energy released in the burst Ei SO with the GRB 
peak energy, E p , and with the break-time of the 
X-afterglow light curve. 

log^o = a+Mog M-t^ +& 2 log - A- (1) 
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where % is measured in days and a and bi, with 
i = 1,2, are calibration constants. 

The other one is the Ghirlanda relation, [6] . It 
connects the collimation-corrected energy, or the 
energy release of a GRB jet E gamma , with the 
peak energy E p 
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where a and b are two calibration constants. The 
Ey term depends directly on the break time tb 
and the isotropic energy in the following way: 
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where Ffc eam — 1 — cos 6*, with 9j et the jet opening 
angle defined by [7]: 
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where E iso _ 52 = E iso /10 52 ergs, n is the circum- 
burst particle density in 1 cm -3 , and r] 7 the ra- 
diative efficiency 

Nevertheless the calibration of these relations 
has been necessary in order to avoid the circu- 
larity problem. This means that all the relations 
need to be calibrated for every set of cosmologi- 
cal parameters. However, currently, there is no a 
sample of low redshift GRBs, up to z = 0.2-0.3, 
to allow a cosmology-independent calibration of 
these relations. In order to overcome this dif- 
ficulty, Liang et al., [8], proposed a method in 
which they calibrated several GRB relations us- 
ing an interpolation method applied to the Su- 
pernova la data. In this way, it becomes possible 
to build a GRB-Hubble Diagram by computing 
the luminosity distance for each GRB by 



di = 



( Ei, 



bole 
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where Ei SO is the isotropic energy emitted in the 
burst and S' bolo is the bolometric fluence corrected 
to the rest frame of the source in consideration. 
The final step consists in calculating the distance 
modulus fi = 25 + 5 log di, and its error, for each 
GRB. 

Our data sample is constituted by 27 GRBs 
taken from the Schaefer sample of 69 GRB, [4]. 
For each GRB we assume the same value for the 
radiative efficiency, ?y 7 = 0.2. Moreover we as- 
sume that the error in the determination of the 
redshift z is negligible, as well as for the same ra- 
diative efficiency. The error bars on the distance 
modulus are photometric uncertainties only, ne- 
glecting any other type of error, e.g. GRB intrin- 
sic variability, so: 



(2.5erio g E iso Y + (1.086a Sbo JS bolo y\ (6) 



with <J\ og E i30 and <7s bolo obtained from the error 
propagation applied to Eq.(fT]) and Eq.Q. We 
note also that the assumptions of a well-known 
no is a strong hypothesis since the goodness of 
the fit depends also on this parameter. For this 
reason we consider the n values for each GRB 
given in the Table 2 of [S], also if it lacks, as 
we said previously, a complete and clear physical 
basis for the considered relations. 
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As a final step we simply add a sample of Su- 
pernova data to our GRB sample, in order to 
constraining better the following data fits at low 
redshifts. We take in consideration the UNION 
catalog, compiled by the Supernova Cosmology 
Project, [10]. 



Table 1 

Results of the fits. LZ is for Liang-Zhang relation, 

GGL for the Ghirlandact al. one 

Fit qa jo + [j go 

LZ -0.68 ±0.30 0.021 ± 1.07 3.39 ±17.13 
GGL -0.78 ±0.20 0.62 ± 0.86 8.32 ± 12.16 



3. DATA FITTING 

So we have obtained a tool capable of testing 
the various cosmological models up to a redshift 
z = 7-8. As a test instrument, we first performed 
an evaluation of the various cosmographic param- 
eters related to Hubble's Series and then we es- 
timated the trend of the equation of state of the 
Dark Energy. The Hubble Series is the Taylor 
expansion of the Hubble law: 



di(z) = d H z\ 1 + - [1 - q ] z 
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We stopped at the fourth order in the expansion 
so that the various terms in the series are related 
with the cosmographic parameters, defined as fol- 
lows: 
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where t represents the cosmic time. These param- 
eters can be expressed in terms of the cosmologi- 
cal density 1— Oa an d the Equation of State (EoS) 
w = p/p of the Dark Energy, [llj . where p and 
p are respectively the pressure and the density 
of the fluid constituting the Dark Energy. Using 
eq.© as theoretical model, we start our data fit 



Table 2 

Cosmological density parameters 



Fit 


Qm 


Q A 


LZ 


0.37 ±0.31 


0.63 ± 1.13 


GGL 


0.28 ±0.30 


0.72 ± 1.09 



in order to estimate the cosmographic parameter 
go, jo and sq. The results of this fit are given in 
Tableffl 

An immediate application of these results con- 
sists in the estimation of the cosmological density 
and Dark Energy EoS parameters. Nevertheless 
we don't know the exact evolution of the Dark 
Energy EoS with the redshift, so that we need 
a parametrization for the function w. For this 
reason we consider the Chevalier-Polarski-Linder 
(CPL) parametrization defined as follows 



w(z) 
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With this function it is possible reconstruct 
the cosmological evolution of the Dark Energy 
and, more important, come back to the epoch 
of the transition acceleration-deceleration of the 
universe. Our previous fit results constrain the 
cosmological parameters to vary in the ranges 
given in Table [21 while for the CPL parameters 
we obtain 



Wq 



^0.53 ±0.64 



0.59 ±0.77, (13) 



that within the errors agree with the ACDM con- 
cordance model but it doesn't agree with the 
epoch of the transition acceleration-deceleration 
: z tr ~ 10. 

This estimate does not agree very well with the 
real value of the Dark Energy EoS, since it is esti- 
mated that the transition redshift is much earlier 
than obtained from our simulations. This is also 
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because the method used here works very well 
only at z < 1. For this reason we need an ana- 
lytical formulation of the Hubble diagram valid, 
in principle, at any redshift, [12] . Let's write the 
Friedmann equation in the case of a quasi-flat uni- 
verse, so with k ~ 0: 



H\z) = Hi (1 + z) 



3(iu+l) 



(14) 



Now the w-parameter indicates the EoS w = p/p, 
where p and p are the pressure and the matter- 
energy density of the Universe, respectively. Us- 
ing the CPL parametrization for w, that in this 
case is the EoS for a fluid representing the total 
Universe, then we rewrite the Friedmann equa- 
tion as follows: 



H{z) = H Q 
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In this way an analytical expression for the lumi- 
nosity distance, and also for the distance modu- 
lus, can be achieved: 



1 + 3u>o + 3w a 3w a 



'2(1 + 



.(16) 



Substituting such an expression in the distance 
modulus, we obtain a model for data fitting which 
could work, in principle, at any z. It is impor- 
tant to stress that the obtained expression for 
the Hubble parameter H(z) is independent of the 
density parameters, Qm and Oa- Moreover it 
is also important to emphasize that we use the 
CPL parameterization not only for the dark en- 
ergy component, as is done in the previous analy- 
sis, but it is used here for the total energy-matter 
density of the Universe, not just because mat- 
ter (dark and baryonic) is contributing by a null 
pressure and the radiation contribute is null, but 
also because we obtain an exact analytical formu- 
lation for the luminosity distance at any redshift 
z. 

With this in mind we start the fit of our com- 
bined sample and the results are given in Table 
[31 Immediately we note that the fit results differ 
if we consider only the Supernova la sample from 



Table 3 

Results of the fits. SNela is just for the Supernova 
la data, LZ is for the GRBs data obtained from 
the Liang-Zhang relation, GGL for the Ghirlanda 
et al. one. Note the improvement on the w- 
parameter using the GRBs data in addition to 
the well-known SNela data, corrected for the 3 
"wrong" GRBs. SNela is just for the Supernova 
la data, LZ is for the GRBs data obtained from 
the Liang-Zhang relation, GGL for the Ghirlanda 
et al. one. 



Relation 


W 


w a 


SNela 


-0.9097 ±0.07 


0.755 ±0.054 


LZ 


-1.39 ±0.38 


1.18 ±0.37 


GGL 


-1.46 ±0.38 


1.36 ±0.32 


LZ + SNela 


-1.15 ±0.10 


0.93 ±0.11 


GGL + SNela 


-1.42 ±0.12 


1.24 ±0.13 



those computed for the sample composed also by 
the GRBs. The residual analysis, see Fig. [21 gives 
us the solution to this problem: a smooth trend 
up to z ~ 3.5 in the data distribution can be 
detected, but beyond this limit we have 3 GRBs 
that exceed the 3tr confidence limit. It seems like 
a sort of darkening at high redshift in the GRB 
luminosity detected. This fact is fundamental for 
the goodness of the fit. Indeed these GRBs rep- 
resent the most distant object in our Hubble Di- 
agram, so their weight on the fit is very high. 
Possible explanations for this anomalous bright- 
ness could be various. In particular there could 
be some process of absorption of gamma radiation 
by low energy photons incoming from the cosmic 
thermal background, [14]. Nevertheless the CPL 
parametrization could be a bad approximation for 
the cosmological EoS so that we should consider 
this problem as due to a uncorrect starting model. 



Anyway we reject these 3 GRBs and repeat the 
analysis without the 3 outliers GRB. The results 
of the fits are given in Table [3] 

From a first analysis of the fit results the com- 
bined, and corrected, GRB+Sn sample is to- 
tally in agree with the concordance model for the 
present epoch. This is confirmed also by a Monte 
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Table 5 

Cosmological density parameters, with uncertainties computed to la confidence limit, obtained from the 
MC-like procedure 



Sample 






n A 


n k 


x 2 


UNION - 


- GRB 


0.26 ±0.14 


0.73 ±0.14 


0.01 ±0.04 


1.032 


UNION - 


- GRB corrected 


0.25 ±0.10 


0.74 ±0.135 


0.01 ±0.035 


1.00027 
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Figure 1. Redshift-Distance modulus diagram for 
the SNeIa±GRB sample. The blue line represents 
the data fit. 
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Table 4 

Results of the fits corrected for the 3 "outlier" 
GRBs. SNela is for the Supernova la data, LZ 
is for the GRBs data obtained from the Liang- 
Zhang relation, GGL for the Ghirlanda et al. one. 

Relation wq w a 

LZ + SNela -0.95 ± 0.01 0.74 ±0.01 
GGL + SNela -0.865 ± 0.005 0.66 ± 0.005 



Carlo- like procedure for the comparison of the re- 
sults with the usual likelihood estimator. In this 
case we consider, for a sake of completeness, also 
a small contribute to the total density of the uni- 
verse due to the curvature density. The results of 
this procedure are given in Table g] and the con- 
tour plot, where the boundaries correspond to lc, 
2cr and 3er confidence levels, are plotted in Fig. 
©. 

We have adopted a similar procedure in the 
case of an EoS evolving with redshift. The re- 
sults of this analysis is plotted in Fig. ([4]) where 
the best fit value, the cross in the figure, cor- 
responds to the value wq = —0.84 ± 0.14 and 
w a = 0.72 ± 0.06, in a good agreement with 
the results obtained using our theoretical model, 
Eq. [16] In particular we obtain for the epoch of 
the transition acceleration-deceleration the value 
z tr = 5.17 ± 0.133, a result that, if higher than 
the redshift of the farther GRB used, could be 
in agreement with current quasar formation sce- 
nario. 



Figure 2. Comparison between the best fit of /z 
and the observed distance modulus /i fcs at any 
redshift. The black dots are the GRBs data and 
the red line is the best fit curve representing the 
theoretical distance modulus. 



4. CONCLUSIONS 

From this analysis we conclude that the cor- 
rected GRB+Sn sample agrees fairly woth the 



6 



— i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — |— 




< ~2, 



- w 1 


| ] L ] L 1 1 1 1 | I I 1 1 | 1 1 1 1 


\ 






¥ 

















O.l 0.2 0.3 



Q.5 0.6 Q.7 



-1.2 -1.1 



-l.O —Q.SM 



-O.T -0_6 



Figure 3. 68%, 95% and 98% constraints on tt m 
and obtained from UNION compilation and 
the GRB sample corrected for the 3 wrong GRBs. 



ACDM concordance model with a small contribu- 
tion in terms of curvature density, being this pa- 
rameter equal to k — 0.01 ± 0.04. In other words 
the theoretical model delineated in Eq. [16] seems 
a good approximation of the observed cosmogra- 
phy and agrees very well with the concordance 
model so that we can argue that GRBs could be 
good distance indicators at redshift values up to 
z = 4. 

However, more robust samples of data are 
needed and more realistic EoS models, with re- 
spect to the simple perfect fluid models, should 
be taken into account in order to suitably track 
redshift at any epoch (see for example p3]). With 
the improving of the observations, in particular 
with the recent launch of new satellites devoted 
to the GRB surveys, as Fermi-GLAST and AG- 
ILE, one should be able to expand the samples 
of GRBs, possibly by data from objects at higher 
redshift. 

In summary, considering these preliminary re- 
sults, it seems that GRBs could be considered 
as a useful tool to remove degeneration and con- 
strain self-consistent cosmological models. Fur- 
thermore the matching with other distance in- 
dicators would improve the consistency of the 
distance-redshift Hubble diagram by extending it 
up to redshift 7 — 8 and over. 



Figure 4. 68%, 95% and 98% constraints on wq 
and w a obtained from UNION compilation and 
the GRB sample corrected for the 3 wrong GRBs. 
The cross represents the best fit value and it is 
in a good agreement with what found using the 
theoretical model described in Sect. 3. 
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